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Abstract 



We study the time evolution of a single spin excitation state in certain linear spin chains, 

as a model for quantum communication. Some years ago it was discovered that when the spin 

rS^ ' chain data (the nearest neighbour interaction strengths and the magnetic field strengths) are 

^pH, related to the Jacobi matrix entries of Krawtchouk polynomials or dual Hahn polynomials, 

so-called perfect state transfer takes place. The extension of these ideas to other types of 

discrete orthogonal polynomials did not lead to new models with perfect state transfer, but 

r^ ' did allow more insight in the general computation of the correlation function. In the present 

O^, paper, we extend the study to discrete orthogonal polynomials of g-hypergeometric type. A 

remarkable result is a new analytic model where perfect state transfer is achieved: this is 

when the spin chain data are related to the Jacobi matrix of g-Krawtchouk polynomials. The 

^ ' other cases studied here (affine g-Krawtchouk polynomials, quantum g-Krawtchouk polynomials, 

OA , dual g-Krawtchouk polynomials, g-Hahn polynomials, dual g-Hahn polynomials and g-Racah 

polynomials) do not give rise to models with perfect state transfer. However, the computation 

of the correlation function itself is quite interesting, leading to advanced g-series manipulations. 



in 

§ '. 1 Introduction 

In quantum information processing, the transfer of a quantum state from one qubit to another 
qubit is a crucial ingredient. For long-distance quantum communication, it is generally accepted 
KN ' that this should rely on optical means. For short-range or mid-range quantum communication 
j^ ■ (e.g. from one quantum processor to another) it seems to be more natural to use qubit chains 
(spin chains of interacting fermions) as quantum wires. In such a chain, qubits interact with their 
nearest neighbours and the interaction is described by some tridiagonal interaction matrix. It 
is advantageous to accomplish state transfer by just letting the system evolve by itself, without 
dynamical control. Bose was the first to propose such spin chains as a quantum channel for quantum 
state transfer between qubits located at different ends of a spin chain tlii2j, see ^ for a review. 
Bose's work inspired many researchers to study the possibility of perfect state transfer in spin 
systems [3HTT]. 

The transmission of quantum states is commonly performed by a chain of qubits coupled via the 
Heisenberg or the XY interactions |12ftl4j . In this context, one often assumes to have individual 
control of the nearest-neighbour couplings in the spin chain, leading to the idea of pre-engineered 
qubit couplings [IKiS] . An interesting consequence of well-chosen controlled couplings is that one 
can obtain mirror inversion of a quantum state with respect to the center of the chain, and this can 
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lead to perfect state transfer [HEllTT] at certain specified times over arbitrary length of the spin 
chain. 

Theoretical investigations of spin chains for perfect state transfer can be divided in two classes: 
analytical solutions and numerical solutions. The original two systems described in [5j are analytical 
solutions: the coupling strength Jk at position k of the system (fc = 0, 1, . . . , A^ — 1) consisting of 
A^+1 interacting qubits is some analytic function of k and N , and there exist closed form expressions 
of the eigenvalues and eigenstates of the single- fermion states. In the case of [5j, these eigenstates 
are related to discrete orthogonal polynomials (namely Krawtchouk polynomials and dual Hahn 
polynomials). Another analytical solution was given by Shi ei a/ [7j. On the other hand, one can use 
a numerical procedure known as the inverse eigenvalue problem, and design spin chains for perfect 
state transfer numerically [9]. In this case, one starts from a (spatially symmetric) set of values 
(the single- fermion energies), and constructs numerically a tridiagonal matrix having these values 
as eigenvalues. The off-diagonal elements of this matrix then yield the values for the strengths J^. 

The analytical solutions have a number of advantages, certainly from the mathematical point 
of view. Apart from the fact that one has a complete analytic description of the single-fermion 
eigenvalues and eigenstates (and thus, through Slater determinant, of all n- fermion eigenstates), 
some other interesting quantities can be computed in closed form. For example, the transition 
amplitude at time t of an excitation from site s to site r (or the "correlation function") can be 
determined analytically [16j. 

In the two elegant analytical solutions of [5j, the single-fermion interaction matrix is related to 
the Jacobi matrix of a system of discrete orthogonal polynomials (Krawtchouk and dual Hahn). 
This inspired Chakrabarti and Van der Jeugt [16] to investigate other systems for which the inter- 
action matrix coincides with the Jacobi matrix of a system of orthogonal polynomials. Although 
no other solutions with perfect state transfer were found this way, the theoretical analysis gave 
rise to a number of interesting results: explicit formulae for transition amplitudes (or correlation 
functions), an explanation of why these two systems discovered in [5] are so special, and a group 
theoretical approach of the problem. 

Apart from several known systems of discrete orthogonal polynomials of hyper geometric type 
in the Askey-scheme [T7] , there is also a list of discrete orthogonal polynomials of g-hypergeometric 
type. In this context, it is a natural question to ask whether Jacobi matrices of these g-orthogonal 
polynomials could also function as interaction matrices for spin chains, and whether they would 
give rise to new solutions with perfect state transfer. This is the topic treated in the current paper. 

Among the main results is indeed a new analytical solution for a spin chain with perfect state 
transfer. This new solution occurs in the context of g-Krawtchouk polynomials. We have studied 
all cases of the q- Askey-scheme, and our analysis has shown that this is the only new case with 
perfect state transfer. Therefore, we shall concentrate on this g-Krawtchouk case, and present all 
details related to this. For the other cases the computations and results are rather complicated, 
and we shall only give a brief summary. 

2 The qubit chain as a spin chain 

In this short section, we shall describe the notation and necessary ingredients for a qubit chain. 
This is a system of A^ -|- 1 interacting qubits (spin 1/2 particles) in a quantum register, with an 
isotropic Hamiltonian of XY type: 

N-l ^ N 

fc=0 fc=0 



where J^ is the coupHng strength between the qubits located at sites k and A; + 1, and /ifc is the 
"Zeeman" energy of a qubit at site k. To describe the Hilbert space associated with the Hamilto- 
nian, one adopts a standard fermionization technique [18j. The Jordan- Wigner transformation [T9] 
maps the Pauli matrices to spinless lattice fermions {a^, a|,| /c = 0, 1, . . . , A^} obeying the common 
anticommutation relations, and in terms of these the Hamiltonian ([1]) takes the form: 
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This describes a set of A^ + 1 fermions on a chain with nearest-neighbour interaction (hopping 
between adjacent sites of the chain), and subject to a non-uniform background magnetic field 
denoted by /i^ (A: = 0,1,..., A^). Initially, the system is in its completely polarized ground 
state |0) = |00---0) = |0) (g) |0) (g) • • • (g) |0), where |0) denotes the spin down state. Let |A;) = 
|00 • • • 010 • • • 0) = a|.|0) (fe = 0, 1, . . . , A^) denote a state in which there is a single fermion at the 
site k and all other sites are empty, i.e. \k) describes the state in which the spin at the site k has been 
flipped to |1). Clearly, the set of states |A;) {k = 0, 1, . . . , A^) forms a basis for the single-fermion 
states of the system. In this single-fermion basis, the Hamiltonian H takes the matrix form 
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(3) 



The dynamics (time evolution) of the system is completely determined by the eigenvalues e-,- and 
eigenvectors ipj of this matrix. It is then a standard technique [5l|T8] to describe the n-fermion 
eigenstates oi H (n < N) using the single-fermion eigenstates ipj and Slater determinants, which 
is why we concentrate on the single-fermion eigenstates. 

The matrix M in ([3]) is real and symmetric, so the spectral theorem [20] implies that it can be 
written as M = UDU^ , where D is a diagonal matrix and U an orthogonal matrix: 



D = diag(eo,ei, 



,e7Vj 



UU^ 
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(4) 
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The entries of D are the single-fermion energy eigenvalues, and the columns of the matrix U are 
the (orthonormal) eigenvectors of Af, i.e. the single-fermion eigenstates ipj = ^,^.=0 ^'^i 1^) '^ith 

The dynamics of the system under consideration is described by the unitary time evolution 
operator U{t) = exp{—itH). Assume that the "state sender" is located at site s of the spin chain, 
and the "state receiver" at site r (s and r are site labels, belonging to {0, 1, ... , A^}). At time t = 
the sender turns the system into the single spin state |,s). After a certain time t, the system evolves 
to the state U{t)\s) which may be expressed as a linear superposition of all the single spin states. 
So the transition amplitude of an excitation from site s to site r of the spin chain is given by the 
time-dependent correlation function 



frAt) = ir\Uit)\s). 
Using the orthogonality of the states (fj , one finds [IB] : 
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One says that there is perfect state transfer at time t from one end of the chain to the other 
end when \fN,o{t)\ = 1- The conditions for perfect state transfer can quite easily be described 
in terms of the "mirror symmetry" of the matrix M in ([3]), see [9l[TT]. However, if our aim is to 
study analytical solutions, we should also require the conditions that the eigenvalues €j and the 
eigenvector components Ukj should be analytic (closed form) expressions. 

We shall now consider the cases where the values characterizing the system (J^ and hk) are 
related to the Jacobi matrix of a set of discrete orthogonal polynomials of g-hypergeometric type. 
This has the advantage that the quantities €j and Ukj are known explicitly, and from these the 
correlation functions fr,s{t) can be computed. Our interest goes beyond perfect state transfer only: 
the aim is to study those cases where the correlation function has a closed form expression. Of 
course, the cases with perfect state transfer deserve extra attention. 

3 g-Krawtchouk polynomials and perfect state transfer 

The purpose of this section is to describe in detail the first case of an interaction matrix of the 
form (l3|) related to the Jacobi matrix of a finite system of discrete orthogonal polynomials of q- 
hypergeometric type. We shall recall some of the necessary notation for g-series. For the case of 
qi-Krawtchouk polynomials, we give the eigenvalues ej and eigenvectors ^j. Then the purpose is 
to compute the correlation function ([7]). Due to the nonlinearity of ej (with respect to j), this 
function turns out to have a complicated structure. Only for certain forms of the deformation 
parameter g, and at specific times t, one can simplify the expression of the correlation function 
fr,s{t)- Fortunately, when g is a positive rational number (quotient of two odd integers), and when 
the other parameter p appearing in the g-Krawtchouk polynomials takes the special form p = q~^ , 
the system yields perfect state transfer. This situation is summarized in subsection 3.3. We end this 
section by reconsidering the general correlation function, and apply some (/-series manipulations in 
order to write it in an appropriate form. 

3.1 Notation for g- functions and g-Krawtchouk polynomials 

The standard reference book on g-hypergeometric functions is |21j . and we follow the notation from 
this book. For a list of orthogonal polynomials of g-hypergeometric type, see |17j . 

In the context of g-series, g is a positive real number with g 7^ 1, and for us it can be considered 
as an extra parameter in the model. We use the common notation for g-numbers: 

N = W, = \^ {n G Z) (8) 

and [n] — )• n in the limit g — )■ 1. For any complex number a and any nonnegative integer n, the 
g-shifted factorial is defined by 

(a; q)n = (1 - a)(l - ag) • • • (1 - ag""!), (9) 

and the product is just 1 when n = 0. Sometimes, in the context when < g < 1, one also uses 

the infinite product 

00 

{a-qU = J{{l-aq^). (10) 

fe=0 

For products of g-shifted factorials, it is common to use the abbreviation 

(oi, 02, . . . , ca; q)n = (ai; g)„(a2; g)n • • • [a-A] q)n- (11) 



The g-hypergeometric series or basic hypergeometric series a(Pb depends on A numerator pa- 
rameters tti, B denominator parameters bi and a variable z and is defined as [21j: 
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In most of tlie series considered here, the numerator contains a parameter of the form q~"^, with 
m a nonnegative integer. In that case, the series (112p is terminating: it has only m + 1 terms 
since {q~^',q)n = for n = m + l,m + 2, . . .. If one of the denominator parameters is of the 
form q~ (with A'^ a positive integer), then one of the numerator parameters should be of the form 
q~"^ with m < A'^ in order to make sure that the series terminates before one reaches zeros in the 
denominator. This will always be the case here. 

Let us now consider g-Krawtchouk polynomials Kn{q~^',p,N;q), characterized by a positive 
integer parameter N and a positive real parameter p: p > 0. This polynomial of degree n in q~^ is 
defined as ^7\ 



Kn{q''')^Kn{q-^-p,N;q)= 3(^2 
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The Q'-Krawtchouk polynomials satisfy a discrete orthogonality relation, namely 

N 
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where the weight function is 



/ N {q ^;q)x, .-x 



and the square norm takes the rather complicated form 
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It is easy to see that (i„ > for < g < 1 and also for g > 1 (since p > 0). The polynomials 
Kn{q~^) also satisfy the following three term recurrence relation: 



(1 - g-")K„((?-^) = AnKn+l{q-'') - {An + Cn)Kn{q-n + A„i^„_i(g-^), 
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It is appropriate to introduce orthonormal g-Krawtchouk functions 
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and then the corresponding orthonormal recurrence relation is 
- [-x]K„(g"^) = -Jn^iKn^iiq-'') + hnKn{q' 



JnKn+liq 
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where 
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Again, it is easy to verify that J„ > (n = 0, 1, . . . , A^ — 1) and /i^ > (n = 0, 1, . . . , N). Following 
the technique of |16j and [22], we now have 

Lemma 1 Let Mqx be the tridiagonal {N + 1) x (N + 1) -matrix (Jacobi matrix) 
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where Jn and hn are given by (|20p . and let U be the {N + 1) x {N + l)-matrix with elements 
Uik = Kj{q-^). Then 
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Herein, 
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In other words, the eigenvectors of the Hamiltonian (in the single- fermion case) corresponding 
to the quantities (|2Up have components equal to normalized (7-Krawtchouk polynomials, and the 
corresponding energy eigenvalues are ej (j = 0, 1, . . . , A^). 

Finally, note that working with all positive J„-values as in ([3]), or with a matrix like (|21|) where 
all off-diagonal elements are negative, does not make an essential difference (see the remark at the 
end of section 1 in [TB] ) . The main difference is that there are also sign changes in the components 
of the corresponding eigenvectors. 

3.2 Computation of the correlation function 

Consider now a spin chain for which the values of J„ and /i„ are fixed by (|2U|) . For this system, the 
eigenvalues of the single-fermion states are given by (p^ . and the eigenvectors by ^j = ^^^=0 ^kj \k) 
with Ujk = Kj{q-''). 

In this spin chain, consider the transition from a sender site s to a receiver site r. The transition 
amplitude or correlation function, given in general by ^, becomes: 
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The purpose is to investigate whether one can compute the sum in ()25p . For this, let us first 
concentrate on the special case when sender and receiver are at different ends of the chain, namely 
when s = and r = N. After some simplifications, one finds for the summation part: 
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Because of the factor z~^~^\ the sum in (|26p is not of g-hypergeonietric type, hence there is no 
hope that it can be simphfied any further for arbitrary values of z (i.e. for arbitrary values of t). 
For certain specific values, however, simplification does take place. Since z = e~**, one has 

Assume now that the deformation parameter (7 is a rational number of the following form: 

p 
q^ =7^) with P and Q odd positive integers (having no common factors). (27) 

Then, for each value of the index k with k < N: 
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Suppose now that we consider the system at time 



t = T = Q^TT, 
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then 

In this case, the expression (126p simplifies drastically, since 
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This is in fact a consequence of the g-binomial theorem {2H (II. 4)]: 
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Taking into account the expressions for dg ^^'i dj\f from ()16p . one has 
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It is not too difficult to see that this expression takes its maximum value when p = q ^ , and in 
that case 
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and thus we have discovered a new analytic model with perfect state transfer. 

Note that at time t = 2T one has z~^~^' = 1 for all k, and thus from (|30p we have JnaC^T) = 0. 
More generally, in that case it simply follows from the orthogonality of g-Krawtchouk polynomials 
and (f25]l that /r,s(2T') = 5r,s- Clearly, the system is periodic in time with period 2T. 



3.3 Perfect state transfer related to g-Krawtchouk polynomials 

It is convenient to collect all the ingredients for the new model in this subparagraph. Perfect state 
transfer takes place if the parameter p appearing in g-Krawtchouk polynomials is equal to q~ . 
p = q~^ . In that case, w{x) and dn take the simpler form 



w{x) 



{Q'-,Q)n x(x-V)l2 



{Q;Q)x{q;q)N-x 
and 

, _^iq,-q;q)n{q,-q;q)N- 



After some straightforward simplifications, one finds for the spin chain data J„ (n = 0, 1, . . . , A^ — 1) 
and hn{n = 0,l,...,N) from (|20]) : 
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and 

Note that /i„ = h^^n (n = 0, 1, . . . , A) and J„ = Jn-i-ti (n = 0, 1, . . . , A — 1), so the matrix (pTj) 
is mirror symmetric, as required for perfect state transfer. 
We can now state the following result. 

Theorem 2 Let q ^ 1 be a positive number of the form q = Q/P with Q and P positive odd 
integers. For the spin chain system ([2]) characterized by the couplings ()3ip and the /i„ 's ()32p . there 
is perfect state transfer from, site to site N at time t = T = Q^tt. At time 2T the system is back 
to its original state (at time 0); in fact the system is periodic in time with period 2T. The single 
fermion eigenvalues are 

ej =q-^ + q'^ + ... + g-^ {j = 0, 1, . . . , iV), 

and the single fermion eigenstates are given by ipj = X]n=o ^nj \n) where 
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3.4 General remarks 

The special case where perfect state transfer is possible occurs when q~^ = P/Q, with P and Q 
odd (positive) integers. What happens in other cases? Clearly, when P and Q would be both even 
integers, they have a common factor that can be canceled. So there remains to see what happens 
when q~^ is a rational number of the form 

even odd 

odd even 

In that case, one can write: 

,-=2'^ or ,-' = 2-'^. (33) 



where r is a positive integer and P and Q are again odd (positive) integers. In the first case, we 
have 
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so the parity of the expression in brackets becomes independent of k. In the second case, one has 



q ^ +q '^ + --- + q '' 



= 2rkQN ^ ^^^ °^^ integer), (35) 

so again the parity of the expression in brackets becomes independent of k. As a consequence, one 
can never obtain that z^^''' = (— l)'^ (or z^'-" of the form x^) for a certain time t in the cases (j33p . 
and this is a condition for further simphfication of the sum (I30p . 

Now that is has become clear that q = Q/P, with Q and P odd integers, plays a special role, 
let us go back to the general case (i.e. arbitrary p) with a parameter q of this form. The expression 
of the general correlation function is given by (j25p : 
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Just as in (j29p . this simplifies at time t = T = Q tt, since then z~^~ ' = (— 1) . The previous 
summation can then be simplified using the (symmetric) Poisson kernel for Al-Salam-Chihara 
polynomials (see |231 Proposition 4] or [24' (14.8)]): 
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Herein, §1^7 is the notation for a very- well-poised s4'7 series {211 § 2.1]. In the current case of 
what results is a terminating §07 series to which Watson's transformation formula [2H (III. 18)] can 
be applied. This transformation formula rewrites the special §07 series as a terminating 4^(j)3 series, 
and after some simplifications one finds: 
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This is a terminating balanced 4(^)3 series, and it can be written in several forms according to Sears' 
transformation [21, (III. 15)]. (The form given here is the most symmetric one, but ifl + A^ — r — s 
is a negative integer smaller than —r or —s, it should be transformed so that the termination of the 
series is not caused by the denominator.) Expression (1381) is a simple formula for the correlation 
function with arbitrary r and s, arbitrary parameter p, but valid only at time t = T (half the 
period). 



Let us now consider again the case that p = q ■ Then (j38p yields 

f frj.. _ 1 (-g~'; q)r{-q~''; g)s(i; q)N-r-siq~^; q)r+s , . 
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So, due to the appearance of (1; q)N-r-s one finds 

frAT) = Sr+s,N. (40) 

In other words, in the case q = Q/P and p = q"^ , there is also perfect state transfer from site s 
to site A^ — s at time t = T = Q tt. 

4 Spin chains related to other orthogonal polynomials of g-hyper- 
geometric type 

The g-Krawtchouk polynomials described in the previous section are not the only q-deforma- 
tions of ordinary Krawtchouk polynomials. In fact, there are four possible (7- generalizations, 
see [17]. The remaining three cases are known as the affine g-Krawtchouk polynomials, the quan- 
tum g-Krawtchouk polynomials and the dual g-Krawtchouk polynomials. All these polynomials 
satisfy a discrete orthogonality relation with finite support {0, 1, . . . , A^}. Apart from these q- 
generalizations of Krawtchouk polynomials, there are three other sets of orthogonal polynomials of 
g-hypergeometric type with a discrete orthogonality over a finite support, namely q'-Hahn polyno- 
mials, dual (7-Hahn polynomials and g-Racah polynomials. In the context of the present paper, all 
these cases should be investigated following the techniques and ideas of section 3. The purpose of 
the current section is to give an overview of the main results for all these cases. First of all, the 
central (but negative) conclusion is that none of these new cases gives rise to an analytic model 
with perfect state transfer. In this sense, the example of section 3 is unique. In view of this, 
we think it is appropriate to describe not all the details of the computations involved for these 
remaining six cases. Although some of the computations of correlation functions lead to interesting 
g-series manipulations, especially for mathematicians fascinated by g-series, we have chosen not to 
overload this section with such details. The reader interested in the actual g-series computations 
should find sufficient information in the summary given in this section. Others can just skip to the 
final expression of the transition amplitude f^fi for transfer from site to site A^ for each of the 
following cases, and establish that perfect state transfer is not possible. 

4.1 AfRne g-Krawtchouk polynomials 

Let A^ be a fixed positive integer, and q a positive real parameter with q ^ 1. The affine q- 
Krawtchouk polynomials K {q~^;p,N;q) are characterized by a positive real parameter p, and 
defined by |17j 






■i(P2 



g-,(7-,0 

N ;?;' 



n = 0, 1,...,A^. (41) 



The affine g-Krawtchouk polynomials satisfy a discrete orthogonality relation, namely 

N 
Y,w{x)K^^{q--)K^%--) = d^6mn, (42) 

x=0 

where [T7] 

= (j^tMt,,)n ,^ ^ (.;cz).(g;g). . ,.^ 

{q;q)x{q;q)N-x {pq;q)n{q;q)N 
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For positivity, one needs 

<p < q~^ when < q <l and < p < q^^ when q > 1, (44) 

and we shall assume that p satisfies this condition. The orthonormal affine g-Krawtchouk functions 



K^iQ-n 



'w[x) 



d„ 



i^r(9-") 



satisfy the following recurrence relation, which follows from [17, (3.16.3)] 



x]^^ (g--) = -J^^,K^^,{q~-) + /ini^r (<?"" 



JuK^^q-n 



with 



«^n 



n ■ 



N]{l-pq 



,n+l^ 



dn 



hr. 



[n]pq 



n-N 



[n-iV](l-pg"+^) 



(45) 



(46) 



(47) 



where J„ > (n = 0, 1, . . . , A^ — 1) and /i„ > (n = 0, 1, . . . , A^). The same result as Lemma[T]now 
applies, but with the spin chain data given by (j47l) . 

Let us now consider the computation of the transition amplitude in a spin chain governed by 
the quantities ([47]) (z = e"**): 



frAt) 



1 



N 



\/drd, 

1 



:^w{k)K^^{q-')Kt%- 



-k] 



** fc=0 

N 



{pq;q)k(.q;q)N{pq) ^ 



Vdrds ^ {q;q)k{q;q)N-k 



■3<P2 



.g-^o 



pq,q 



-N iq^q 



392 



q ^,q *,0 
pq,q~^,0 



-k] 



(48) 



As in subsection 3.2, this sum cannot be simplified unless q~^ is of the special form q~^ = P/Q 
with P and Q positive odd integers, and t = T = Q^ir. For this time t, we have z~[~*^l = (—1)'^, 
and now one can continue with the computation of ()48p . In this case, the product of the two s4'2 
functions is computed by applying the product formula of g-Hahn polynomials given by [211 (8.3.3)] 
and substituting a = p and 6 = in this expression. Following this, the product of the two 302 
series in (f^Sl) can be rewritten as a sum (X^m=o) ^^^^ terminating 31^3 series (terminating because 
q~"^ appears as a numerator parameter in the 3i;^3). Exchanging summations over k and over m, 
and performing the inner summation over k (which can be simplified using the g-binomial theorem) 
finally leads to the following: 



frAT) 



(-i;g) 

\/drd. 



N . 

N_sr^ {q 

^ m=0 



'-'',q'-'';q)m{pq'^T"' 



{q,-q 



l-N 



-N.. 



■3</'3 



N-r+l-m N-s+l~m ^q^Pq 

pq^ q J q 



2N-m+3 



(49) 

The current expression can in general not be simplified further, except in the special cases when 
the sender is at one end of the chain (s = 0). One finds: 



and in particular 



/r,0(T) 

fNAT) 



(-1;' 



)N 



Vdrdo 

-i;g)Af 



2<?'l 



q^-^,0 1 

—q pq^ 



{-i;q)N{pqf'^Vimq)^- 



(50) 



(51) 



It can be verified that this expression never reaches the value 1 (for A^ > 1), due to the condi- 
tions (I44p for p. So perfect state transfer is not possible in this case. For certain values of p and g, 
so called high-fidelity transfer can still be achieved |lj, and in this context such models might be 
worth considering. 
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4.2 Quantum g-Krawtchouk polynomials 

The quantum (7-Krawtchouk polynomials K'i^"^^{q~^;p,N;q) are characterized by a positive real 
parameter p, and defined by [17] 



K^'^'iq-n = K''^{q''';p,N;q) = 2^1 



-N 1I1PI 



The orthogonality relation is of the form ()42|) , but now with 

{pq; q)N~ 



n = 0,l,...,A^. 



(52) 



i(;(x) 



_^_^Y^x{x~l)/2^ 



{q;q)xiq;q)N- 

j _ (g; Q)N-n[q,Pq] q)n , ^ \N-n^N ^Nn+N{N+t)l2~n(n+\)l2 

"n — / X \ ^) p q 

[q,q;q)N 

For positivity, one needs 

p > q~ when < q < I and p > q~ when g > 1. 



(53) 



(54) 



The orthonormal quantum gr-Krawtchouk functions satisfy the following recurrence relation, which 
follows from [17, (3.14.3)] 



where 



[-x]K^'^iq-n = -Jn-iKf^M'^ + KKt'^iq-'') - .InK^^Aq-^, (55) 

[n-N] fd^ ^ [n]{l-pq^) [n - N] 



Jn 



hr. 



-,2n 



-,2n+l ■ 



(56) 



pq2n+i Y d^ ' pq^"- pq^ 

The computation of the transition amplitude in a spin chain governed by (I56|) is again a rather 
technical question. As in the previous subsection, the summation part cannot be simplified unless 
z~^~^' = (—1)'^, and this happens when q~^ = P/Q with P and Q positive odd integers, and 
t = T = Q^TT. Using similar techniques as in subsection 4.1, one finds: 



N 



frJT) = (-1) 



N{-'^;q)N {pq;q)N s^ {q '',q ''■,q)m{pq 



E 



r+s+l—N\m 



Vd^s {q;q)N ^^ {q,- 



/^-^,g-^;g)m 



„r—N „s—N „—m 

1 '9 '9 .„ „-l„JV-m+2 

„-l„-Ar r-m+1 s-m+l ^q^P </ 



X 303 

The current expression can in general not be simplified further, but one has 



fN,,{T) = (-1; q)N p-''q-^''''^''^/'^{-ir{pq; q)N. 



(57) 



(58) 



Due to the conditions (I54p this expression never assumes the value 1 (for A^ > 1), so perfect state 
transfer is not possible. 

4.3 Dual g-Krawtchouk polynomials 

The dual g-Krawtchouk polynomials K{\{x)] c, A^; q) are polynomials of degree n in \{x) = q~^ + 
(,qX-N ^ characterized by a real parameter c with c < 0, and defined by [IT] 



Kn{\{x)) = Kn{\{x)-c,N-q) = sh 



,q ■^,cq 



x-N 



;q,q 



n = 0,l,...,N. (59) 
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The orthogonality relation is of the form (j42p , but now with 



w[x) 



(eg ^,q ^;g)^. 1-cg^^ ^ _^ xi2N-x) , _ (g;g)n(c ^ g)jv , _jvNn .«nA 



The orthonormal dual g-Krawtchouk functions satisfy the following recurrence relation, which fol- 
lows from dZl (3.17.3)] 

- [-x]{l - cg"-^)i^„(A(x)) = - J„_iK„_i(A(x)) + hnKn{X{x)) - JnKn+i{X{x)), (61) 

where 

Jn = -[n- N]^l^, K = -N - cq-^'in - N]. (62) 

So a similar result holds as in lemma [H except that the single particle eigenvalues are now given 

by 

efc = -[-k]{l - cq''-^) = (g-i + g-2 + . . . + g-fc)(l _ cq''-^). 

Now one can compute the transition amplitude in a spin chain governed by (j62p . As in the previous 
subsections, the summation part cannot be simplified unless 2;^[~"(^~^'' ' = (— 1)'^. This happens 
when q~^ = P/Q with P and Q positive odd integers, and moreover c = —2Q (or an integer 
multiple of this), and for the time t = T = Q^tt. Just as in subsection 4.1, one can use (a limit of) 
the product formula for g-Hahn polynomials in order to rewrite the product of dual g-Krawtchouk 
polynomials as a sum over certain 3i;^2-series. Next, one has to change the order of summation. For 
the inner sum over k there appears, apart from various g-shifted factorials, the quotient ~^ _j^ 
(due to the weight function), and this has to be rewritten as 

l-cg^"-^ _ (cgi-^;g)2fc _ (cV2g^-^/^ -cV^g^-^/^; g)^ 
1-cg-^ (cg-^;g)2fc {cy^q-^/\-cy\-^/^;q)k ' 

Then, the inner sum reduces to a special case of a very-well-poised g^s-series, which can be summed 
according to [21^ (11.20)]. The final result is: 

^"^'^' ^ - Vd;d; (cgi-^; g2)^ ^^ (g, -gi-^v, q-N, ,)^ ^^^ '^ M c) 



-m{m-l)/2+{N-r-s)m. ^ , 



N-r-m+l N-s-m+1 ' "/' '-^H 

y 5 q 



(63) 



The case (r, s) = {N, 0) yields the simplest expression: 



(_l;g) ^(_c)^/2g-7V (JV-l)/4 

(cgi-^;g2)^ 






Again, one can verify that this expression never assumes the value 1 (for A^ > 1), so perfect state 
transfer is not possible. 

4.4 q-Hahn, dual g-Hahn and g-Racah polynomials 

The remaining g-hypergeometric discrete orthogonal polynomials (with a finite support) are more 
complicated, and the notation becomes quite heavy. We shall give some details regarding the g- 
Racah polynomials, and then some of the results for g-Hahn or dual g-Hahn polynomials will follow 
by taking an appropriate limit. 
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Let A^ be fixed and assume < g < 1. The (7-Racah polynomials Rn{n{x);a,l3,'y,5]q) are 
polynomials of degree n in fi{x) = q^^ + j5q^^^, and involve four parameters a, /3, 7 and 5: 



Rn{fi{x)) = Rn{fJ.ix);a,(3,'j,6;q) = 4(^)3 



q-'',al3q''+\q-=',j6q''+^ _ 

aq,l36q,-fq ' ' 



(65) 



where n = 0,1,..., A^. One of the following relations should hold: aq = q~^ , fiSq = q~^ or 
75 = q~ . Here, we shall assume that pSq = q~ , so in the following we put 5 = /3~^q^ . 
The g-Racah polynomials satisfy the orthogonality relation: 



N 



'^w{x)Rm{Kx);a,(3,^,6;q)Rn{fJ.{x);a,l3,-/,6;q) = dnSmn, (66) 



x=0 



where 



and 



{aq,(35q,jq,j6q;q)r, 1 - -^Sq + 

w{x) = _i ■- ^_i T ^ z ^ [af3q) , (67 

{q,a ^-fdq,p ^7q,0q;q)x I - -fSq 

^ {a/3q'^,Pj-^;q)N {q, a/37~^g» a/3g^+^ /3g; g)^ 1 - «/3g .^,1 „Ar.„ , . 

" (a/?7-ig,/3g;g)jv ((?-^,ag,a/3g,7g;g)„ 1 - a/3g2n+i ^^ ^« ^- ^ ^ 

Certain conditions should hold for positivity of the weight function, e.g. < a < q~^, < /? < q^^ 
and 7 > g~^ (then automatically 6 > q~^); another option is a > q~^ , (i > q^^ and < 7 < 
q~^. The orthonormal g-Racah functions Rn{fJ-{x)) = Rn{fJ-{x))y^w{x)/dn satisfy the following 
recurrence relation (see [T71 (3.2.3)]) 

- [-X](l - 7<5g^+l)i?„(/i(x)) = -Jn-lRn-lif^ix)) + KRniKx)) " JnRn+liK^)), (69) 

where 

Jn = -r^\[^, hn = -^^^, (70) 

I- q\ dn I- q 

with 

_ (1 - ag"+i)(l - a/3g"+i)(l - /3(5g"+i)(l - 7g"+i) 

"~ (l-a/3g2n+l)(l_a/3g2n+2) ' 

_ g(l - g")(l - /3g")(7 - a/3g")(^ - ag") , . 

(l-a/3g2")(l-a/3g2"+i) " ^' ^ 

We have once again a similar result as in lemma [H but with the single particle eigenvalues given 

by 

The next item to work out is the computation of the transition amplitude in a spin chain gov- 
erned by (j70p . As in the previous subsections, the summation part cannot be simplified unless 

2/3 = (— 1) . This can occur when q = P/Q with P and Q positive odd integers, 

and moreover 7 = 2Q^ P (or an integer multiple of this), and for the time t = T = Q^ir. The 
actual computation of fr,s{T) then involves the following steps: 

• In the sum 

N N 

^wik)Rsifiik))Rrifiik)) = Y,Hk)Rk{Ks))Rk{Ks)) 

k=0 k=0 

the product of two g-Racah polynomials Rk{fJ-{s))Rk{fJ-{r)) is written as a single sum over 
io(/'9 series, using f2T] (8.3.1)]. 
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-[-k]il - j5q'+^) = (g-1 + g-2 + . . . + g-'=)(l - Igfc-^) 



• Exchanging the order of the two summations, the inner sum can be performed using the very- 
well-poised 6</'5 summation theorem |21[ (11.21)]. This finally leads to a single sum expression 
of io(j)Q series. 

This last expression reads: 



[P IQ ,-i',q)N sr^ [q ,q ,a p q ,a p q 



f rrp-s ^ \p m ,-^,q)N sr^ 



) y;m 



Vd;d:{/3-^^q^-^-q^)N ^^ {q,a-^p-^^q-^ , l3-^q-^ ,a-^l3-^q-^-\-q 



l-N a-N.„\ 



xg IP iq ,q)m.u^, ^afiqN--+^ , - ^a/Jg^— +1 ,aq,^q, a(3q''+'' , a/3g^+'^+2— , 



a/3^-ig^-'^+i, g-™, q-r^ q~s ^ aPq''+^,a(3q'+^ 

g N+s+2-m N-r+l-m N-s+l-m 



l^,„N+2 



,f3-'iq 



(72) 



This is quite an impressive expression. Note that when the sender is at site {s = 0), the iQt/'g 
series collapses, and the remaining sum reduces to a terminating 4^3 series. In the simple case that 
(r, s) = (A^, 0), this simplifies further and yields: 



f ,rj.. ^ (-l;g)^(7//?)^/^g-^(^-^)/" l i-l)^{aq,f3q,jq,af3^-^q;qU .„.. 

^''''^ ' W-'7q'-'';q^)N V («/3<?2,a/3g^+i;g);v ' ^ ' 

Unfortunately, one can again verify that this expression never reaches the value 1, so perfect state 
transfer is not possible. 

Observe that the g-Hahn polynomials Qn{q~^]Oi,l3;q) are obtained from the g-Racah polyno- 
mials in the limit 7 — )• 0. Without giving all details here, let us mention that in this case the 
spectrum is of the form (|24l) . and that for i = T = Q'^tt (where q~^ = P/Q, P and Q positive odd 
integers), the correlation function is obtained from ()72p by taking the same limit. In particular, 
one has in this case: 

In a similar way, the dual g-Hahn polynomials i?n(^(x); 7, 5, A^; q) are obtained from the g-Racah 
polynomials in the limit a — )• (and using /3 = q~^~^/"f)- The spectrum is of the form 

efc = -[fc](l-7<5(z'+'). 

Again under special conditions {q~^ = P/Q, P and Q positive odd integers; 7^ = 2P , t = T = 
Q^tt), the correlation function is obtained from ()72p under the same limit. For the special case 
{r,s) = {N,0), one finds: 

InAT) = /ja'^l'y J{iq,Sq;q)Nhqr. (75) 

A particular simple expression is obtained in this case if both of the parameters 7 and 6 of the dual 
g-Hahn polynomials are equal. Then ([75]) yields 

fN,o{T)^ =l^^{7qr^'. (76) 

However, for the allowed range of parameters (for 0<g<l: 0<7< q^^ and < 5 < q^^; or 
7 > q~ and 5 > q^ ) the expression is always less than 1. So once again, perfect state transfer 
is not achievable. 
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5 Summary and conclusions 

In this paper, we have considered linear spin chains with a nearest-neighbour hopping interaction 
as models for quantum communication. We have studied the time evolution of single fermion 
states in such a spin chain. Certain special spin chains allow perfect state transfer [HEIIH], and 
these systems can be related to Krawtchouk polynomials or dual Hahn polynomials. Following 
the success of these two initial systems, it was a logical step to study spin chain models related 
to other discrete orthogonal systems. This was in fact the topic of an earlier paper [16j. In that 
paper, no new models for perfect state transfer were discovered, but the orthogonal polynomial 
approach allowed the explicit computation of the correlation function (or transition amplitude) for 
the known models. Furthermore, the models were also approached from a group theoretical point 
of view. 

The next logical step is then to study spin chain models related to discrete orthogonal polyno- 
mials of g-hypergeometric type. This was exactly the topic of the present paper. The main feature 
is that the spin chain data Jt and hk in ([2]) are now given not by ordinary numbers but by certain 
(/-numbers (of course, since g is a positive parameter, we end up with real numbers again once q 
is fixed). We have now investigated all possible spin chain models with an interaction matrix ^ 
coinciding with the Jacobi matrix of a discrete orthogonal system of g-hypergeometric type. The 
fascinating outcome is that one new spin chain model for perfect state transfer has been discovered 
this way, in relation to g-Krawtchouk polynomials. The details of this model and the actual com- 
putations have been described in detail in section 3. A general feature of all the models given in 
this paper is that the single fermion eigenvalues ej are no longer linear in j (as for the Krawtchouk 
case, see 0]) or quadratic in j (as for the dual Hahn case, see [Ij). Instead, the single fermion 
eigenvalues are of the form ej = q^^ + g~^ -!-••• + q~-' (or a multiple of this). As a consequence the 
correlation function fr,sit), involving summations over e~**"^^ (see dZ])), can no longer be simplified 
in general. However, when 5 is a certain rational number, we have still managed to compute and 
simplify the correlation function for specific values of the time t. It is in this context that the new 
model for perfect state transfer was encountered. 

For the remaining models based on a Jacobi matrix of affine g-Krawtchouk polynomials, quan- 
tum g-Krawtchouk polynomials, dual g-Krawtchouk polynomials, g-Hahn polynomials, dual g-Hahn 
polynomials or g-Racah polynomials, perfect state transfer is not achieved. For all these additional 
models, we have computed the correlation function (in the cases where simplification takes place). 
These computations on their own have some interest, as some remarkable g-series manipulations 
can be performed. Furthermore, these models might still be worth considering in cases where 
high-fidelity transfer is of importance. 

It is worth mentioning the difference with some recent work on quantum state transfer in a 
g-deformed chain [25]. In [25], the algebraic relations for the fermion operators a^, a[, in ([2]) 
themselves are g-deformed. Then, it is rather the representation theory of g-deformed algebras 
that plays a role. Here, the operators in ([2]) are ordinary undeformed fermion operators. But the 
fermion chain data, i.e. the strengths Jk and hk, are given by g-numbers. 
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